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Abstract 

We study at a quantitative level the impact of the uncertainties on the value of the 
W boson mass measured at hadron colliders due to: i) the proton parton distribution 
functions (PDFs), ii) the value of the strong coupling constant a s and Hi) the value of the 
charm mass used in the PDF determination. The value of the W boson mass is extracted, 
by means of a template fit technique, from the lepton-pair transverse mass distribution 
measured in the charged current Drell-Yan process. We study the determination of m w 
at the Tevatron and at the LHC with 7 and 14 TeV of center-of-mass energy in a realistic 
experimental setup. The analysis has been done at the Born level using the event generator 
HORACE and at NLO-QCD using the event generators DYNNLO and ResBos. We 
; consider the three global PDF sets, CTEQ6.6, MSTW2008 and NNPDF2.1. We estimate 

that the total PDF uncertainty on m w is below 10 MeV both at the Tevatron and at 
the LHC for all energies and final states. We conclude that PDF uncertainties do not 
challenge a measurement of the W boson mass at the level of 10 MeV accuracy. 
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1 Introduction 

The measurement of the W boson mass represents a very important test of the Standard Model 
and of its extensions, like e.g. the MSSM, and provides indirect bounds on the mass of the Higgs 
boson p~H3]. This measurement has reached a very high level of accuracy: the current world 
average is m w = 80.398 ±0.023 GeV [I] and the best single experiment measurements have been 
obtained by DO [5j and CDF 0[7] at the Fermilab Tevatron with m w = 80.401 ± 0.043 GeV 
and m w = 80.413 ± 0.048 GeV respectively. The prospects for the combined measurements at 
the end of the Tevatron run, with 4 fb _1 of total collected luminosity, are of a final error of 
roughly 15 MeV [8]. The prospects for the measurement at the CERN LHC are at the level 
of 15 MeV, or even 10 MeV [9j[T0]. At this level of accuracy it becomes necessary to quantify 
in detail the various sources of theoretical uncertainties that contribute to the final systematic 
error. 

The mass of the W boson is measured at hadron colliders in the charged current (CC) Drell- 
Yan (DY) process by studying the charged lepton transverse momentum p\ distribution, the 
missing transverse momentum p v t distribution or the lepton pair transverse mass distribution, 
defined as 

MY = ^^(l-cos^-^)) , (1) 

where the neutrino four-momentum p% and angle <p v are inferred from the transverse momen- 
tum imbalance in the event. The mass of the W boson is obtained by fitting the experimental 
distributions with the corresponding theoretical predictions, where m w is kept as a free param- 
eter. 

A measurement of m w at the 10 MeV level is not only a very ambitious goal from the 
experimental side, but it is also very challenging from the theoretical point of view due to 
the careful modelling of the production mechanism that is required. We can illustrate these 
difficulties with the following example. It is known that the result of a fit of m w to a given 
theory template is very sensitive to the shape of the distributions. In Fig. [T] we consider two 
transverse mass distributions at the Born level obtained with two values of m w which differ by 
10 MeV. If one takes the ratio bin by bin of the histograms, one sees that a small shift of 10 
MeV in m w induces a non trivial distortion of the shape at the permille level. Therefore, if we 
aim at measuring m w at the 10-20 MeV level, we should, from the theoretical side, have the 
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Figure 1: Left plot: the transverse mass distributions at the Tevatron at Born level obtained with 
two values of m w which differ by 10 MeV. Right plot: the bin by bin ratio of these two distributions. 
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Table 1: Cross sections within acceptance cuts, at Born level, as a function of m w . We also show the 
percentage difference between pairs of cross sections that differ by 10 MeV. 



control on all the perturbative and non-perturbative corrections which can change the shape 
of the relevant kinematic distributions at this level of precision. 

On the other hand, the total integrated cross-section is not significantly affected by changing 
m w . As shown in Table (H a shift by 10 MeV of m w yields a change of the cross section at 
the 0.04% level. Thus, it is important to disentangle the normalization effects, which are very 
weakly related to the precise value of m w , from the effects that modify instead the shape of 
the distributions, which have a larger impact on the measurement of m w . 

The Drell-Yan cross-section is given by the convolution of the parton distribution functions 
(PDFs) of the two incoming hadrons with the partonic cross-section. The crucial role of QCD 
corrections to the partonic processes has been widely discussed in the literature [THUS]- The 
very important role of the 0{ot) EW corrections in the precision study of the CC DY process 
is also well known (for a complete list of references see [13]). It is the aim of the present paper 
is to study three different sources of uncertainty related to the PDFs and their impact on the 
measurement of m w : 

1. the PDFs are affected by uncertainties due to the error of the experimental data from 
which they are extracted, as well as by theoretical uncertainties like the non-perturbative 
functional form parametrization. These uncertainties affect the prediction of the DY 
observables and, in turn, the extraction of the value of m w . Moreover in some cases 
the central values obtained from different PDF sets differ more than the nominal PDF 
uncertainties: we need to account for this by considering more than one PDF set; 

2. the NLO-QCD corrections sizably modify the Born level lepton transverse momentum 
distribution and, more moderately, also the transverse mass distribution. The precise 
effect of these corrections depends on the value of the strong coupling constant, which is 



ultimately correlated with the PDFs and with their evolution. Therefore also the precise 
value of a s should be taken into account in a precision determination of m w ; 

3. the PDFs depend on the value of the heavy quark masses m c and m;, due to two different 
reasons: the first one is the fact that 0(m 2 c /Q 2 ) terms have a non- negligible impact on 
PDF fits, and the second that heavy quark PDFs are obtained by assuming them to vanish 
at threshold, and then to be generated by perturbative evolution. For these reasons, 
the value of the m c used in the PDF determination has an impact on the kinematic 
distributions from which m w is extracted, and tuhs must be accounted for. The value of 
nib on the other hand does not affect W production due to the smallness of b-initiated 
contributions. 

Uncertainties related to PDFs are known to be an important component of the total sys- 
tematic error in the determination of m w at hadron colliders. In the most recent CDF and DO 
measurements PDF uncertainties are estimated to be between 10 and 13 MeV [3]. Ref. [TO] 
estimates PDF errors in m w prior to LHC data to be ~25 MeV, decreasing at the few MeV 
level once the constraints from LHC processes are taken into account. On the other hand, 
there are claims [U] that with the current knowledge of PDFs a determination of m w with 
a precision Am w < 10 MeV is far from being possible. In this paper we want to revisit the 
impact of PDFs and related uncertainties on the determination of m w at the Tevatron and the 
LHC, considering the most updated global PDF sets and related theoretical uncertainties, like 
the values of a s and m c . 

The paper is organized as follows. In Sect. [2] we present the general strategy that we will 
follow to estimate the shifts of the measured value of m w induced by PDF uncertainties. In 
Sect. [3] we present the numerical results of our analysis for the transverse mass distribution and 
in Sect. H]the results for the PDF impact on the determination of m w . In Sect. [5] we explore 
the improvements on PDF uncertainties for the determination of the W mass provided by LHC 
data and in Sect. Owe draw our conclusions. 

2 The determination of m w : general strategy 

In this section we present the general strategy that we adopt to estimate the impact of PDF 
uncertainties in the determination of m w at hadron colliders. First of all, we introduce the 
fitting procedure and its validation. Then we discuss the event generators and settings adopted 
to compute the theoretical distributions. Finally, we discuss the PDF sets that are considered 
in this study together with related sources of theoretical uncertainty. 

2.1 The fitting procedure and its validation 

We consider in the present study differential distributions in Charged Current Drell-Yan pro- 
duction generated with different PDF sets and we treat them as samples of pseudodata. The 
Montecarlo error on each bin is taken in the statistical analysis as the error affecting the pseu- 
dodata. The pseudodata are generated with a given common nominal value of m w called m^,, 
which is taken to be = 80.398 GeV, the current world average. 

The general fitting strategy is summarized in Fig. [2j First of all, we generate the templates 
for a given fixed PDF set, in this case the central set of CTEQ6.6, and for different values of 
m w , with very high statistics, IB events at Born level. Then for each member of the PDF sets 
considered, including the error PDF sets, we generate pseudo-data with fixed = 80.398 GeV 
using exactly the same event generator as for the templates, with lower statistics, 100M events 
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Figure 2: Flowchart that summarizes the procedure used to determine the shift in m w induced by 
any given PDF set. More details are provided in the text. 

at Born level. Then we compute the y 2 between the pseudo-data and each of the templates: 
the template with best y 2 provides the information on which is the shift in m w induced by 
this particular PDF set. As expected for consistency when pseudo data is generated with the 
central CTEQ6.6 set, we get y 2 ~ 1 for the selected template with m w = 

The templates have been computed for 100 (at Born level) and 20 (at the NLO-QCD 
level) different values of the W mass. The range for these templates has been taken to be 
80.398 ±0.050 MeV at Born level and 80.398 ±0.036 at NLO-QCD. We compare each template 
with the pseudodata and compute the reduced y 2 function, defined as 

i N hias f/nj /^data\ 2 
2 _ 1 \ ^ \yi u i ) ■ _ -I jvr (n\ 

X-3 ~ AT f-.da.taN 2 ~~ ' ' ' ' ' iV templatcs \*) 
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where Oi is the value of the i— th bin of the distribution O (e.g. the W transverse mass) and 
the superscript refers to the pseudodata or to the j— th template. The value of m w used in 
the template which minimizes y 2 is considered as the preferred value of m w and the difference 
Am w = m w — m^y, is the shift induced by the PDF set chosen for that set of pseudodata. A 
similar approach has been used in [T0~l[T5l [T6"] . 

The fitting procedure has been validated by using samples of pseudodata that have been 
produced with the same inputs and the same event generator of the templates but with different 



2 





/ 1 M 




/ 10 M 




/ 100 M 




/ 340 M 




nominal MW=80.398 GeV 


templates with IB events 


Tevatron 


pp -> W ± -* fiv^ 





80.38 80.385 80.39 80.395 80.4 80.405 80.41 80.415 80.42 

w (GeV) 

Figure 3: x 2 distributions obtained fitting Born level pseudodata with Born level templates for 
the same fixed PDF set at the Tevatron kinematics. The different curves correspond to different 
pseudodata samples each with different statistics. The Ax 2 = 1 rule indicates the resolution, at 68% 
C.L., on the W mass. 

statistics. In this case the function \ 2 denned in Eq. [2] can be used to make a x 2— test. When 
fitting pseudodata obtained with a given event generator, the nominal value of m w used in 
the generation of the data is rediscovered as preferred value within an interval determined by 
the condition A\ 2 = 1, which can be interpreted as a 68% CL. interval. This interval shrinks 
as the number of events considered increases and correspondingly their statistical fluctuations 
are damped, as shown in Fig. [3J In this example the templates have been generated with IB 
events while the pseudodata has been generated with increasing statistics from 1M to 340M 
events. We also checked that, fitting 1000 independent samples of Born level pseudodata, the 
corresponding minima follow the \ 2 distribution, as expected. 

When the statistics of the pseudo-data become close to those of the templates, the \ 2 can 
deteriorate since it becomes sensitive to the statistical fluctuations of the latter, not accounted 
for in Eq. [2J This imposes a practical limit on how accurate the pseudodata can be. This effect 
can be seen in Fig. [3] for the case of pseudo-data generated with 340M events. We find that 
a good compromise between resolution and stability with respect to fluctuation is provided by 
using templates of IB events with pseudo-data generated with 100M events. 

2.2 Event generation 

Let us discuss now how the theoretical predictions of the DY kinematic distributions have been 
generated. We have studied the production process pp — > fi + + X at the Tevatron Run II (s/s 
= 1.96 TeV). We also consider the two processes pp — >■ ji + + X and pp — > pT + X at the LHC 
for -y/i = 7 TeV and \/s = 14 TeV center-of-mass energies. In the absence of QED effects, not 
considered here, our results will be identical to those obtained with electrons instead of muons. 
The numerical results have been obtained using the following values for the input parameters: 



Tevatron 


LHC 


p" ± > 25 GeV 
$ T > 25 GeV 

kl < i-o 


> 25 GeV 
$r > 25 GeV 
M < 2-5 



Table 2: Selection criteria for — > l^v events for the Tevatron and the LHC. 



G„ = 1.16637 10- 5 GeV" 2 m w = 80.398 GeV m z = 91.1876 GeV 

T w = 2.141 GeV sin 2 9 W = 1 - m 2 w /m 2 z m H = 120 GeV 

V cd = 0.222 K s = 0.975 V c6 = 

Kd = 0.975 K s = 0.222 Kb = 

V td = V ts = Vfb = 1 

The charm quark in the partonic cross section is treated as a massless particle, while the bottom 
quark does not contribute because of the vanishing top density in the proton. 

In the generation of Drell-Yan charged current events we used the selection criteria sum- 
marized in Table EJ These kinematic cuts are similar to those used in the corresponding 
experimental analysis. Note that the main difference between the Tevatron and LHC cuts is 
a wider acceptance for the rapidity of the leptons in the latter case. The W transverse mass 
distribution has been studied in the interval 50 GeV< MY < 100 GeV, with a bin size of 0.5 
GeV, since the jacobian peak region is the most sensitive for the determination of m w . All the 
following analysis are performed with bare leptons both in the pseudodata and in the templates. 

The pseudodata and the templates have been generated using the following event generators: 
at Born level with HORACE [Tf], at NLO-QCD with DYNNLO [T8j and at NLO+NNLL-QCD 
with ResBos [19] . These generators allow to compute the distributions of the final state leptons 
in the DY processes at various perturbative orders. For example ResBos includes, on top of the 
NLO-QCD corrections, part of the NNLO-QCD terms matched with the resummation of the 
large log j at leading logarithmic (LL) and next-to-leading logarithmic (NLL) accuracy, 

and has been widely used at the Tevatron. 

Our final results for the determination of m w will be those obtained at NLO-QCD with 
DYNNLO, although, as we will show below, the qualitative results are already very similar at 
Born level. 



2.3 PDF uncertainties 



The proton PDF sets considered in this study are the three global sets that include all the rele- 
vant hard scattering data. In particular we will use the NLO-QCD CTEQ6.6 [20], MSTW2008 [2T] 
and NNPDF2.1 [22] PDF sets. Each collaboration provides a prescription to estimate the PDF 
uncertainties: in particular we recall the formula for the symmetric error in the Hessian ap- 
proach (CTEQ,MSTW) 
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and the average over the ensemble of PDF replicas (NNPDF) 



\ 



N 



(3) 



r£ p k=l 



Nrep \ ' 

iJ2(nti {k) }}-(nti}})) 2 ) ■ 

k=l J 

We refer to the original publications as well as to the recent reviews [23TT25] for more details. 
Let us recall that the use of the three global PDF sets is the basis of the current PDF4LHC 
recommendation [26] for the use of PDFs in the analysis of LHC data. 

On top of the PDF uncertainties that arise from the experimental uncertainties of the 
data used in their determination, there are other sources of theoretical uncertainties closely 
related to PDFs. In the first place, PDFs are correlated with the value of the strong coupling 
constant a s (mz) used in the PDF determination, expecially the gluon PDF. Again, all three 
groups provide prescriptions on how to combine the PDF and strong coupling uncertainties in 
a consistent way. A summary of the prescriptions recommended by each group can be found 
in the PDF4LHC working group interim report [25] (see also Ref. [27]). A practical guide on 
the way to efficiently implement the recommendations by the different groups can be found 
in Ref. |28j. While the impact of variations on the value of a s (Mz) are known to be small 
for vector boson production.], they may need to be taken into account at the level of precision 
required for the determination of m w . 

On top of the value of the strong coupling, PDFs depend as well on the value of the heavy 
quark masses m c and m?, due to two different reasons. The first one is the fact that even 
though most LHC perturbative computation are done up to power-suppressed terms, terms of 
O (m 2 h /Q 2 ) do still have a non-negligible impact on PDF fits, expecially to the HERA collider 
data. Power suppressed terms are accounted for in the various General-Mass VFN schemes 
used in modern PDF sets [301432] . and the choice of m c affects the GM-VFN predictions and 
thus the fitted PDFs. Different GM-VFN schemes have been compared in the Les Houches 
heavy quark benchmark study |33J, elucidating their differences and similarities. 

The second reason has simply to do with the fact that heavy quark PDFs are obtained by 
assuming them to vanish at threshold, and then to be generated by perturbative evolution. 
But changing the mass also changes the position of the threshold, and thus the heavy quark 
PDFs (and their contribution to the cross section) depend on the value of rrih- For example, 
for W production, the initial state with one charm and one strange quarks occurs at the Born 
level approximately in the 7% of the cases at the Tevatron, in the 16% at LHC 7 TeV for W + 
production and in the 25% for W~ production, in the 24% at LHC 14 TeV for W + production 
and in the 32% for W~ production. 

For these reasons the precise value of the m c has an impact on the kinematic distributions 
from which m w is extracted and must be accounted for, expecially since charm mass variations 
are known to induce sizeable effects for W production at colliders [2"2"II3"4] . 




3 PDF uncertainties for the transverse mass distribution 

Now that the setup of the analysis has been presented, we consider how PDFs and related 
uncertainties affect the lepton pair kinematic distributions, in particular the transverse mass 
distribution, and in the next section we will consider their impact on the determination of m w . 
As discussed in Sect. [2], only those sources of uncertainties that induce distortions on the shape 
of the distribution (rather than on its normalization) will have an impact for the extraction of 
m w . 



*As opposed to other relevant LHC processes, like Higgs boson production via gluon fusion, where a s 
uncertainties can be the dominant theoretical uncertainty [291 . 



The transverse mass distribution has the advantage, with respect to the lepton pr distri- 
bution, that QCD-NLO corrections are rather moderate and in particular have a small effect 
on the shape of the distribution. Experimental issues in its measurement like the systematic 
uncertainties due to the neutrino px reconstruction are not addressed here. 

In the following we will consider two related distributions: the transverse mass distribution, 

° ^ = dW ^ ' Mf = \/2^r(l-cos(0<-^)) , (4) 

and the same distribution but normalized to the integrated cross section in the region used for 
the m w fit, 

with MY Mq = 50 GeV and Mf ' min = 100 GeV. The motivation to define 6 is that in this 
way normalization effects, irrelevant for the m w determination, cancel out, and one is left only 
with the contribution of PDF uncertainties that induce shape distorting effects. The use of 
normalized distributions has also been adopted in the Tevatron analysis [I]. 

The same NLO PDFs are used both to generate the Born and NLO-QCD distributions. In 
Figs. H]and|5]we compare, for the three PDF sets, the relative size of the pure PDF uncertainties, 
at the Tevatron and at the LHC 7 and 14 TeV, for the transverse mass distributions computed 
at Born level with the HORACE generator. In the latter case we consider separately the two 
cases of W + and W~ production, since in a proton-proton collider the two distributions are 
different unlike in a proton-antiproton collider. We show both the standard, Eq. HJ and the 
normalized, Eq. [5j distributions. 

We observe that the PDF uncertainties in the normalized distributions are much smaller 
than in the standard transverse mass distributions: the reason for this is that variations in the 
normalization of the distribution, which are not relevant for the determination of m w , cancel 
out in the normalized distributions. Note that from Figs. H] and [5] we see that PDF uncertainties 
are at the few permille level. 

The previous plots show that PDF uncertainties are similar for the three global PDF sets. 
However, it could still be the case that the distributions obtained with the central set of each 
PDF set differ sizeably among them, leading to an uncertainty in m w much larger than the 
nominal PDF uncertainty of a single set. To check that this is not the case, in Figs. |6] (for 
the LHC) and [7J (for the Tevatron) we show the ratio of transverse mass distributions for each 
central PDF set normalized to the central CTEQ6.6 predictions. 

The results of Figs. E] and [7J show that, while the standard transverse mass distributions 
differ at the few percent level between different PDF sets, the normalized distributions on the 
other hand are much more similar, providing an excellent agreement of the central values and 
differing only at the permille level. This is the same order of magnitude as the intrinsic PDF 
uncertainties. This suggests that the determinations of m w from the three different sets are 
consistent within the respective PDF uncertainties: we will explicitly verify this expectation in 
Sect. U 

The uncertainties on the transverse mass distributions, this time computed with the DYNNLO 
generator at NLO-QCD, are shown in Figs. [8] and |9j The QCD corrections introduce a new 
partonic sub-process (qg — > qlv\) and the related gluon density uncertainty. The latter induces 
an increase of PDF uncertainties in the large tail of the transverse mass distribution above the 
jacobian peak, where the cross section steeply falls, as well as for small transverse masses. On 




Figure 4: Relative PDF uncertainties in the Born level transverse mass distributions, computed with 
respect the respective central PDF set. From top to bottom: LHC 7 TeV W + and W~ and LHC 
14 TeV W + and W~ . Both the PDF uncertainties on the standard distribution, Eq. [43 and the 
normalized distribution, Eq. are shown. 

the other hand, in the region near the peak, most relevant for the determination of m w , the 
PDF uncertainties at NLO-QCD are similar to those of the Born distributions. 

In the transverse mass distributions normalized to their respective cross sections, the differ- 
ence in PDF normalization has been removed and the uncertainty is due only to the different 
shapes induced by the PDF sets considered. The comparison in Figs. [8] and [9] shows that the 
typical size of the PDF uncertainty on these normalized observables is well below the 1% level, 
whereas in the non- normalized case it ranges between 2 and 3%. The latter are the typical 
PDF uncertainties that one finds for the inclusive cross section [22J. 

In Fig. [10] we show, in the case of the NNPDF2.1 set, how the PDF uncertainties in the 
NLO-QCD transverse mass distribution varies with the energy, collider type and final state. 
The different uncertainties are very similar in size (e.g. they are all at 2% level, below 80 
GeV). Fig. [TU] shows that PDF uncertainties in the transverse mass distribution are relatively 
independent of the collider and final state. This result is reassuring since it shows that, at least 
from the PDF point of view, the determination of m w at the LHC is not more challenging than 
at the Tevatron. 

Let us now assess the impact of the uncertainties related to the values of a s and m c on 
the transverse mass distribution. In Fig. [TTJwe show for NNPDF2.1 the PDF-only uncertainty 
compared to the combined PDF+a s uncertainty. Following Ref. [27] we assume that the un- 
certainty on the strong couping is 5 as = 0.0012 at the 68% confidence level. For simplicity we 
show only the distributions at the LHC 7 TeV: the distributions for Tevatron and LHC 14 TeV 




Figure 5: Same as Fig. |3]for the Tevatron. 



are quantitatively very similar. We conclude that a s uncertainties are negligigle as compared 
to the PDF uncertainties for this distribution. 

We have also studied the dependence of the results on the value of m c used in the PDF 
determination (using the NNPDF2.1 set with m c variations), taking fully into account all 
correlations between m c and the PDFs. In Fig. [T2] we show the ratio of transverse mass 
distributions computed with different m c in the PDFs, divided by the results of the central 
NNPDF2.1 set. It is clear from these results that a different choice of the charm mass in the 
evolution of the parton densities yields a different overall normalization of the transverse mass 
distribution, but it affects very moderately the shape. This is confirmed by the normalized 
distributions: the percentage difference with respect to the reference m c value is consistent 
with zero within statistical fluctuations. 

In summary, we found in this section that PDF uncertainties in the transverse mass distri- 
bution can be kept at the permille level by normalizing them to the integrated cross section 
in the fitted interval. These PDF uncertainties turn out to be very similar for all colliders, 
energies and final states, and are in reasonable agreement between different PDF sets. The 
theoretical uncertainties related to m c and a s , that are important for inclusive cross sections, 
turn out to be negligible for the normalized kinematical distributions. In the next section we 
will assess the impact of these various uncertainties on the determination of m w . 



4 PDF uncertainties in the determination of m w 

We have shown in the previous section how PDF uncertainties distort the shape of the transverse 
mass distribution. Now we use the fit setup presented in Sect. [2] to extract for each template 
the associated value of m w , and check how the values of m w obtained with different PDF 
sets compare to each other and with their intrinsic PDF uncertainties. We fit the W mass 
separately with each different Montecarlo replica (for NNPDF2.1) or with the various Hessian 
eigenvectors (for MSTW08 and CTEQ6.6), and then apply the corresponding prescriptions to 
compute the best estimate for m w and the associated PDF uncertainty for each set. 

In Table [3] we present the results obtained when fitting the Born level normalized transverse 
mass distributions (Eq. \5§. Then in Table @] we show the analogous results obtained when 
fitting the standard distributions (Eq. H]). In both cases we quote the intrinsic PDF error from 
each set, denoted by 5 p df (in GeV), as well as the shift between each set and the reference value 
obtained with CTEQ6.6, denoted by A p df (in MeV). 
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Figure 6: Relative difference between the distributions obtained with the central PDF set of CTEQ6.6, 
MSTW08 and NNPDF2.1, normalized to the CTEQ6.6 result. We show the results both for the 
normalized and for the standard distributions for LHC 7 TeV (upper plots) and 14 TeV (lower plots). 



We note that the values of m w obtained with the standard distributions are shown here 
only for illustration of the sensitivity of the template fit procedure to the normalization choice. 
The templates have been prepared at Born level, separately for each energy, collider type and 
final state. In this way we can claim that in each case we are probing only the effect due to the 
PDF uncertainty. We remark that the central value of CTEQ6.6 coincides, by construction, 
with the value {m^ = 80.398 GeV) used when generating the pseudodata. 

Let us consider first the results obtained with normalized distributions, shown in Table El 
The central values obtained with MSTW2008 and with NNPDF2.1 (that is the spread of A pdf 
values) differ at most by 6 MeV with respect to and lie in general in a A p df ~ 2 — 4 MeV 
interval. The PDF uncertainties are stable when considering different colliders, energies and 
final states. If we now look at the results obtained with the standard distributions, Table HI we 
observe that the central values are spread in a larger interval (±10 MeV) and that also the PDF 
uncertainties are correspondingly increased, <5 p df ~ 5 — 8 MeV. However, it is remarkable that 
even for the standard transverse mass distributions PDF uncertainties turn out to be rather 
small and similar for all colliders and final states. 

In Table [5] we present the results obtained by fitting the transverse mass distribution gener- 
ated at NLO-QCD with DYNNLO, with different PDF sets. The main difference with respect 
to the study at Born level comes from the gluon contribution, absent in lowest order. We fit 
the distributions, normalized to their cross-section in the fitting interval, using DYNNLO tem- 
plates prepared with the central CTEQ6.6 set and normalized. In Table [5] we also provide the 
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Figure 7: Same as Fig. [6] for the Tevatron. 
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Table 3: Results for the determination of m w from normalized transverse mass Born distributions. 
We show in each case the central value of the fit of m w and the spread due to PDF uncertainties, 
<5 p df in GeV. We also indicate well as A p( jf (in MeV), the shift in central predictions from each set 
compared to the CTEQ6.6 reference. 



average (x 2 ) obtained from the fit to the error PDFs in each case. These results are represented 
graphically in Fig. [13j 

The estimate of the PDF uncertainties for m w is quite stable at different energies and 
colliders. The values obtained here are moderately larger than at Born level. To estimate the 
total PDF error, following the PDF4LHC recommendation, we take the envelope of the results 
from the three different PDFs sets. This total PDF error obtained with the envelope method, 
denoted by 5*°^ in Tabled is always smaller than 10 MeV (see also Fig. [TBI . Note in particular 
the excellent agreement of the results from the three sets at the Tevatron, both in central value 
and in PDF uncertainty, yielding a total PDF uncertainty of only 5^ d \ = 4 MeV. 

The above estimates of the PDF uncertainties are somewhat smaller than previous esti- 
mates: for example, Ref. [10] estimates 5^ ~ 25 MeV prior to LHC data@ We would like to 
emphasize that the key in reducing the PDF uncertainties is fitting to the normalized kinematic 
distributions, in a way that normalization effects in PDF uncertainties, irrelevant for the m w , 
cancel out. To illustrate this point, we note that from the results for m w obtained from the 
Born standard distributions, Table SJ and using the envelope of the three PDF sets, one finds 
81% = 12 MeV at Tevatron, <5g = 7(6) MeV at W + {W~) LHC 7 TeV and <5J& = 12(9) MeV at 
W + (W~) LHC 14 TeV, larger than the results of Table [3] (5 MeV, 6 MeV, 5 MeV, 5 MeV and 
4 MeV respectively) and closer to previous estimates. Our estimates for PDF uncertainties on 
m w at the Tevatron are also smaller than existing CDF and DO estimates jl]. There is work in 



2 This uncertainty, computed from CTEQ61, is to be understood as a 90% C.L. and thus ~ 15 MeV at 
68 % C.L. 




Figure 8: Same as Fig. U]but now the transverse mass distributions have been computed at NLO-QCD 
using the event generator DYNNLO. 

progress trying to understand these various results and the differences and similarities of the 
approaches. 

In Table [6] we present the results obtained by fitting the normalized transverse mass dis- 
tributions, at the Tevatron, with ResBos and with CTEQ6.6. The templates used in the fit 
have been computed with ResBos as well, with the central CTEQ6.6 set. Note that within the 
public version of ResBos only the CTEQ6.6 set can be used. By construction the central values 
of the fit coincide with the nominal input value m^. The results for the PDF uncertainties are 
similar to those obtained with DYNNLO at NLO-QCD. 

Let us consider now the impact of the values of m c and a s on the determination of m w , 
which we know to be small from the analysis of the transverse mass distributions. In Table [7] 
we show the results found in the determination of m w for the case of the NNPDF2.1 fits 
with varying strong coupling. As expected from the distributions of Fig. [TTJ, differences are 
negligible, confirming that the uncertainty in a s does not play a role for the determination of 
m w at hadronic colliders. 

In Table [8] we show the results obtained by fitting the Born level transverse mass distribu- 
tions generated with the NNPDF2.1 sets extracted with different values of m c . As it is evident 
from Fig. [T^l the results differ mainly in their normalization. Indeed we observe that the results 
of the Born level fit of the normalized distributions have very small deviations with respect to 
the reference value for the charm mass, at most 4 MeV. This is to be compared with the sizeable 
shifts, at the percent level, observed in the inclusive cross sections when m c is varied in the 
global fit. Again the origin of these different behaviors is that while cross sections depend on 
the normalization of the distribution and thus of the PDFs, m w depends only on its shape. 




Figure 9: Same as Fig. [5]but now the transverse mass distributions have been computed at NLO-QCD 
using the event generator DYNNLO. 
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Figure 10: The relative PDF uncertainty in the standard transverse mass distributions for NNPDF2.1 
for different colliders, energy, and final states. 

5 The impact of LHC data on the m w measurement 

In the previous sections we have discussed the impact of PDF uncertainties on the determination 
of m w at hadron colliders. The PDF sets considered there summarize our understanding of 
the proton structure prior to the LHC. However, LHC data is already providing stringent 
constraints on available PDF, and thus it will further reduce the effects of PDF errors in the 
extracted value for m w . As an illustration of this point, in this section we explore the impact 
that recent ATLAS and CMS measurements of the lepton asymmetry from W decays have on 
the determination of m w . 

The CMS and ATLAS experiment have recently presented their measurements of the lepton 
charge asymmetry from W bosons at 7 TeV. The ATLAS analysis [35] corresponds to muon 
asymmetries while the CMS analysis contains both the electron and the muon asymmetries 
for two cuts of the lepton transverse momentum, px > 25 GeV and pt > 30 GeV. We have 
determined the impact of these LHC measurements on the NNPDF2.1 set by means of the 
Bayesian reweighting technique of Ref . [37] . Preliminary results were presented in [38] , and a 
more detailed analysis will be presented elsewhere. For the comparison with LHC data, the 
theoretical predictions were generated at NLO with the DYNNLO event generator with the 




Figure 1 1 : Comparison of the PDF-only uncertainty and the combined PDF+a s uncertainty of the 
transverse mass distribution for NNPDF2.1. For simplicity we show only the distributions at the LHC 
7 TeV, the distributions for Tevatron and LHC 14 TeV are quantitatively very similar. 
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Table 4: Same as Table [3] for m w fits to the standard transverse mass distributions. 



same binning and cuts as in the respective experimental analysis. 

We consider two cases: one in which the NNPDF2.1 set includes the impact of the published 
CMS and ATLAS lepton asymmetry data and another in which the NNPDF2.1 set includes the 
impact of hypothetical future measurements of the same asymmetries with a relative accuracy 
of 1% (the average error of the published data is about 7%). 

We show the m w distributions obtained with the N rep = 100 replicas of NNPDF2.1 at the 
LHC 7 TeV and the same NNPDF2.1 replicas reweighted with the lepton asymmetry data 
in Fig. [T3] (left plot). The spread of the distribution indicates the PDF uncertainty in the 
determination of m w . We can see that present data already act in the direction of narrowing 
the distribution thus reducing the PDF uncertainties in m w , although the constraints are still 
moderate. 

Larger effects are expected in the scenario with lepton asymmetry pseudo-data with a 1% 
bin-per-bin total experimental uncertainty. We can see that these very accurate pseudo-data 
have the potential to narrow the m w distribution and thus to reduce the PDF error on m w 
by a factor of two or even more. PDF uncertainties could be further decreased if additional 
observables, sensitive to the quark and antiquark combinations relevant for m w production, 
were considered. One example is provided by the accurate measurement of the Z boson rapidity 
distribution at the LHC, that would constrain the small-x sea quarks. 

This exercise confirms that, though PDF uncertainties in the determination of m w are 
already small, they can be further decreased systematically by LHC measurements. 




Figure 12: For NNPDF2.1 we show the dependence on m c of the transverse mass distribution, 
expressed as relative deviation from the central NNPDF2.1 set with m^. =2 GeV 2 . We show results 
both for the normalized and for the standard transverse mass distributions. We consider only the 
LHC 14 TeV case, where charm mass effects are known to be more important. 
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Table 5: Results for the determination of m w from normalized transverse mass NLO-QCD distribu- 
tions. We show in each case the central value of the fit of m w and the spread due to PDF uncertainties, 
<5 p df in GeV. In the right column of each PDF set, the average {x 2 ) P er degree of freedom obtained in 
the fit of the PDF error sets is shown, as a measure of the fit quality. For each collider and final state, 
the final column estimates the total PDF uncertainty 5^ using the envelope method, as discussed in 
the text. A graphical representation of the results is shown in Fig. [T3j 



6 Conclusions 

In this paper we have presented a detailed study of the impact of PDF uncertainties on the 
accurate determination of the W boson mass in hadronic collisions. We have concentrated on 
the shape of the transverse mass distribution and we have used a template fit technique to 
determine a preferred m w value, isolating the PDF effects from other sources of theoretical 
uncertainties. 

Our main conclusions are the following: 

• The Born level study shows that the prediction of the central values and of the PDF un- 
certainties agree between the different PDF sets and are stable when comparing different 
colliders, energies and final states. 

• The NLO-QCD study shows results analogous to the Born level case, with a moderate 
increase of the PDF uncertainty induced by the gluon initiated subprocesses. 

• The use of accurate templates, prepared for each specific collider, energy and final state, 
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Figure 13: Graphical representation of the results of Table [5] for the various colliders and final states 
considered. In each case we draw the envelope of the results from the PDF sets to define the total 
PDF uncertainty (±<5*°f f ) as a thick solid line. The dashed line marks the position of the nominal 
value m^, = 80.398 GeV used to generate the pseudo-data. 
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Table 6: Results for the determination of m w from normalized transverse mass NLO+NLL QCD 
distributions generated with ResBos. We show the central value of the fit of m w and the spread due 
to PDF uncertainties, £ p df in GeV for the case of the Tevatron. The distributions, normalized to the 
corresponding cross section in the fitting interval, have been computed using ResBos and have been 
fit with templates prepared as well with ResBos (with the central set of CTEQ6.6). 



allows to disentangle the role of the PDFs, while keeping fixed all the other input param- 
eters. 

• A sensible and more accurate fit of the W mass can be obtained by studying the shape of 
kinematical distributions, removing normalization effects which should not be explained 
in terms of m w shifts. 

• PDFs and related uncertainties (a s ,m c ) are estimated to be smaller than 10 MeV at the 
LHC for all energies and final states, even before accounting from the improvements from 
LHC data. This implies that PDF uncertainties will be smaller than other systematic 
uncertainties. 

• PDF uncertainties, that are already rather moderate, can be further reduced using LHC 
data alone, without the need of a new dedicated experimental program to constrain PDFs. 
We have illustrated this point using the recent lepton asymmetry data from CMS and 
ATLAS. Measurements of the Z rapidity distribution and other observables will soon 
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Table 7: Central value of the fit of m w obtained with NNPDF2.1, using PDF sets that differ by 
the a s (mz) value, for different colliders and energies. The fit has been done on normalized distribu- 
tions and using normalized templates, and the distributions have been generated at NLO-QCD with 
DYNNLO. 
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Table 8: Central value of the fit of m w obtained with NNPDF2.1 sets with different values of m c for 
different colliders and energies. We include the default value in NNPDF2.1, m? c = 2 GeV 2 as well. 



further reduce PDF uncertainties. Therefore a measurement at the level of 10 MeV 
precision at the LHC, while challenging from many other points of view, does not seem 
to be forbidden by the uncertainties in our knowledge of the proton structure. 

The precision determination of m w is one of the goals of the current 7 TeV run at the LHC, 
due to its potential to indirectly probe new physics at the electroweak scale. This study ensures 
that an accuracy of 10 MeV is certainly within reach, at least in what concerns our present 
knowledge of the structure of the proton. 
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